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Abstract. We present a new exact solution for the twist of an asymmetric thin 
elastic rods. The shape of such rods is described by the static Kirchhoff equations. 
In the case of constant curvatire and torsion the twist of the asymmetric rod 
represents a soliton lattice. 
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The study of elastic rods QIEIU] has a long history and is a subject of increased 
interest 0IE1E1IZ1IH1 especially in connection with bio-mathematical models of proteins 
and of DNA [HlEIl- Thin elastic rods are described by the set of Kirchhoff equations 

On the other hand the static Kirchhoff rod model is also related to the dynamcs 
of spinning tops. In the prersent letter we will consider only the statics of thin elastic 
rods which is described by the static Kirchhoff equations. 

Here we show that when the curvature and the torsion are constant, the twist of 
a thin elastic asymmetric rod (i.e. with asymmetric cross section) can be nonlinear 
along the rod. This represents a new sector of integrable/ exactly solvable cases for 
the Kirchhoff rod. 

Until now it has been assumed that the twist of an asymmetric elastic rod is 
always constant [H]. Now it is clear that the curvature of the rod plays an essential 
role and it may "interact" with the asymmetry of the cross section in order to produce 
a new class of highly nontrivial (soliton lattice) solutions. 

The Kirchhoff rod represents in general a non-integrable system [21 01 E]- 
Nevertheless there are special cases which are integrable: where the twist along the 
rod and the torsion along the cental line of the rod are constant. The case when the 
torsion and the twist are constant along the rod is widely discussed in the literature 

A thin rod can be modelled by a space curve R(s) joining the loci of the centroids 
of the cross sections together with the (so-called) Frenet frame (t(s), n(s), b(s)) 
attached to the rod material and evolving along the shape according to the Frenet- 
Serret equations (with an arc length parameter s) 

dt dn / \ , , \ dh 

— = Kn(s), — = -Kt(s) + rb(s), — = -Tn(s) (1) 
as as as 

Here ti{s), ni{s) and bi(s) are the tangent, the normal, and the binormal vectors 
to the curve respectively. The tangent vector is a unit vector given by t = (^), 
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the curvature k{s) of the curve at the point s is then given by: k{s) := |^| and 
the normal and the binormal vectors can be determined from The local basis 
(di(s), d2(s), d2(s)) which is attached to the rod can be expressed through the Frenet 
frame as follows: 

di(s) = — n(s) sin + b(s) cos (2) 
d2(s) = n(s) cos 4> + b(s) sin 0, ^si^) — t(s) 

where (p is the twist angle of the rod. The components of the derivatives of the local 
basis (d3(s), d2(s), di(s)) with respect to s can be expressed by using the twist vector 

k(s) = Kldi + K2d2 + Ksdg: 

^ =k(s) X d,(s), z = 1,2,3. 
as 

The static Kirchhoff equations describe the shape of the rod under the effects 
of internal elastic stresses and boundary constraints, in the absence of external force 
fields. If F(s) denotes the tension and M(s) denotes the torque of the rod, then in 
the generic case the torque is 

M(s) = Ki(s)di(s) + aK2(s)d2(s) + 6K3(s)d3(s) (3) 

The constant a measures the asymmetry of the cross section and b is the scaled 
torsional stiffness. In particular for symmetric {a — 1) hyperelastic {b — 1) rods 
we have M(s) = k(s). The elastic energy of the Kirchhoff rod is given by: 

H=- M kds=- (Kl+anl + bul) ds. (4) 

The conservation of the linear and angular momenta is provided by the static Kirchhoff 
equations: 

^ = 0, ^ + d3(.)xF(s)^0. (5) 
as as 

The twist vector k is related to the curvature /t(s) and to the torsion r(s) of the center 
line of the rod by: 

k(s) = (fci, fc2, fcs) = (Ksin(/>, KCOS0, r + (f>s), (6) 

where the twist angle is a function of the arc length parameter s: (j) = 4'is)- Then 
the vector equations Q, projected onto the local basis (di,d2,d3) give the following 
6 scalar equations for the tension F = (i^i, F2, F3) and the twist k = (^1,^2,^3): 

+ K2F3 - = (7) 

as 
dFo 

+ Kai^i - >iiFz - (8) 

as 
dF-i 

+ niF2 - K2F1 = (9) 

as 

Fi ^ + (b ~ l)KiK3 (10) 

as 

i^2-^ + (6-a)^2«3 (11) 
as 

6^ + (a-l)Ki/«2-0. (12) 
ds 

In this system of equations we replace fci, ^2 and from eqn. ©. 
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In the case of constant curvature k{s) = fco ^-nd torsion r(s) = tq the equation 
p2|l reduces to the famous static (scalar) sine-Gordon equation for the twist angle: 

^ + ^^iT^^o sin^/(s) = 0, uis) = 2^(s). (13) 

This second order differential equation represents in itself a completely integrable 
Hamiltonian system and allows so-called "soliton"-like solutions. It appears in a wide 
variety of physical problems for e.g. charge-density-wave materials, splay waves in 
membranes, magnetic flux in Josephson lines, torsion coupled pendula, propagation 
of crystal dislocations, Bloch wall motion in magnetic crystals, two-dimensional 
elementary particle models in the quantum field theory, etc. 

The nontrivial solutions of eq. H13(l are single solitons or kinks and soliton lattices. 
We will show that the kink-solutions of (|13() are not compatible with the full set of 
Kirchhoff equations. 

In what follows we will prove that the periodic (soliton-lattice) solution of (|13|l 
is compatible with the full system Q - (|12() and thus represents an exact solution for 
the thin rod in the case of constant curvature and torsion. This solution is given by 
the following expression }12\ : 



(s) — 2 arccos 




(14) 



where m is the modulus of the Jacobian elliptic function sn, and K{m) is the complete 
elliptic integral of the first kind. In the limit m — > 1 K{m) oo and the period tends 
to infinity. The period of the solution is 
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—AmK[m). 

1 — a fco 

In order to prove the compatibility of H13|l with the full system of Kirchhoff equations 
l|7jl- l|12(l we follow the ideas of JJj. Let us consider the following linear combination 
of IZj) and © taking into account the explicit parametrization (0) of the twist vector 
k: 

[Fus - {to + (t>s)F2 + fco-Fa cos </)] sin </> + [Fa^s + (ro + 0,,)^2 - fco-Fa sin 4>] cos (j) = 0(15) 
which reduces to 

Fi^s sine/) -I- F2^s cos4>+ (tq -I- (t)s)[Fi cos(f>- F2 sin(/)] = 0. (16) 

Here and below the subscript s means a derivative with respect to s. Substituting 
l(TT|l and (O into (HHJl yields 

[46-|-2(l-a)cos2<?!)]<?!)^, = (1 - a)(20, + tq)^ sin2(/), (17) 

for the case of nonvanishing curvature fc. Then we use the expression for (j)ss from 
H12|l and write the above equation in the following form 

(1 - a)fc^ sin 40 + 2b[2kl - Acjil - 4to0s - t^] sin 20 = 0. (18) 
After multiplication of both sides of H12|l by ips and integration we get: 

b(j)l ^ {I ~ a)kl sin^ + Co, (19) 
where Cq is an integration constant. Now we eliminate 0^ from (|18|l and (|19|) and get 

4bTo(f>s ^ 3(1 - a)k^ cos 20 + 2{b + a- l)fc^ - 4Co - br^. (20) 
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Next we square both sides of the last equation and ehminate (j)^ using eq. H19|l . Then 
eqn. (|20|l transforms into a pohnomial for cos 20: 

A cos^ 2(p 



C = 0, 



(21) 



bT, 



br'^]'^ 



(22) 

OJ, (23) 
(24) 



B cos 2(j) 
The coefficients take the form 

A = -9(l-a)2fc^, 

B = - 86to2(1 - a)kl - 6(1 - a)kl[2{b + a - l)kl - ACq 
C = 86r2[(l - a)kl + 2Co] - [2(6 + a - l)^^ - 4Co 

If we replace in eqn. (|21|l cos 20 with sin 20 and cos0 with sin0 we will cast the 
trigonometric polynomial in the following form: 

- As\T?2(t>~2Biim^ (t> + A + B + C ^G. (25) 

The above expression must vanish identically for 0(s) from eq. (|14() if this solution is 
to be compatible with the full set of Kirchhoff equations. First let us note that from 
eqn. I)14|l satisfies the following relation jE] 

^"''■^ ■ ' Co, (26) 



" ; 2 • 2 

-kn sm 



where Co = "^-jp^o ('m") ■ ^^'^ '^^^ ^q. (|26|) to replace sin^ and sin^ 20 in eqn. 
H25|l . The expression in H25|) vanishes identically if the following two conditions are 
fulfilled: 

B = -2.4, A + B + C . (27) 

\rn J 

Inserting the first into the second equation leads to: 

2" 

/ m' \ 

1 



C^A 



3(^ 

m 



(28) 



These two conditions are easily fulfilled. The first equation leads to a condition on the 
integration constant Co and the second one with the use of the first gives an algebraic 
equation of fourth order for the curvature fco Solving ean. (|27(l for Co gives: 



Co 



m 
m 





'1-a 6" 


/ 12 


[ 4 +2J 



«,g. 



(29) 



We note here that Co must be positive for the nontrivial case ko ^ and to 7^ 0. 
Finally the condition H28|) leads to the following equation: 



27 



(6-a)(26-a+l)A:;^ + - 



46(26 - a+l) + ^6^(1 -a) 



ki - -6^x4 = 0. 



This equation has always positive solutions for fcp and thus we have shown that 0(s) 
from (|14|l is compatible with the full set of static Kirchhoff equations. 

As we noted before eqn. Ijl^fl has also a kink solution 0(s) = arctan e"/^" with 
So — ■ For this solution it is easy to show that the constant Co in H26() must 

be zero. This is only possible ii ko — and tq = (see eq. H29(l ^. In this case however 
the equation H12|l turns into a linear ODE for 0(s): 0ss = 0. 

Now we turn our attention to the elastic energy density for the thin rod. Wc will 
show that the new exact nonlinear solution for the thin elastic rod may be energetically 
more favorable than the trivial no twist solution = 0. 



An Exactly Solvable Case for a Thin Elastic Rod 



5 



The energy density (per unit length) of the rod is: 

h=^[b{To + (j)sf + {l~a)klsin'^(t) + ak^] (30) 

The Eulcr-Lagrangc equation for the energy is once again the eqn.(13) and its 
nontrivial solution is given in eqn.(14) with sq = The corresponding static 

energy per soliton of the soliton lattice is given by: 

Stwist = l{bT^ + akl) + 2ro[<^(0 - <t>{-l)] + (^("^) ' ^(m')'^(m)) , (31) 

where E{m) is the complete elliptic integral of second kind and m' is the 
complementary modulus: m' = Vl — m^. 

Now let us compare this energy .Btwist with the energy of the trivial solution: If 
there is no twist {(j) = 0, tt, 27r, ...) the energy for length I which corresponds to one 
period d of the soliton lattice 

-E„o twist = Kbro + akl) (32) 

For some choice of the parameters ko,To, a, b, -Etwist may turn smaller than -Eno twist- 
This may happen if the following condition holds true (depending on the choice of tq, 
kq, a, b and m): 

2rom) - <^(-0] + (^("^^ - ^("^')^i^(m)) < 4nro - l^^^d < 0, 

where d is the period of the soliton lattice. We have used that |(/)(/) — (p{—l)\ < 2it. 
In this case the static energy per unit length of the solution lattice will be lower than 
the static energy in the trivial no twist case. 

With this new exactly solvable case for the thin elastic rod the following general 
picture for the twist of thin rods emerges: 

(A) The twist for all symmetric thin rods (a = 1) with constant torsion tq = const 
is either constant (p{s) = const or periodic and linear in s: (f){s) = const. s. Here the 
twist is disconnected from the curvature. 

(B) The twist for all asymmetric thin rods (a ^ 1) (with constant curvature 
Kq and torsion tq) is either constant: (f) = mr/2 or periodic and nonlinear in s: 
(j){s) — 2 arccos[sn(.s/so, ?Ti)]. Here the asymmetry (1 — a) couples the twist with the 
curvature kg which leads to this nonlinear solution. 
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